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A central problem in any quantum theory of gravity is to explain the emergence of the 
classical spacetime geometry in some limit of a more fundamental, microscopic description of 
nature. The gauge/gravity-correspondence provides a framework in which this problem can, 
in principle, be addressed. This is a holographic correspondence which relates a supergravity 
theory in five-dimensional Anti-deSitter space to a strongly coupled superconformal gauge 
theory on its 4-dimensional flat Minkowski boundary. In particular, the classical geometry 
should therefore emerge from some quantum state of the dual gauge theory. Here we 
confirm this by showing how the classical metric emerges from a canonical state in the dual 
gauge theory. In particular, we obtain approximations to the Sasaki-Einstein metric underlying 
the supergravity geometry, in terms of an explicit integral formula involving the canonical 
quantum state in question. In the special case of toric quiver gauge theories we show that our 
results can be computationally simplified through a process of tropicalization. 
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to explain the emergence of the classical spacetime geometry 

in some limit of a more fundamental, microscopic description 
of Nature. The AdS/CFT-correspondence (or “gauge/gravity 
correspondence”) introduced in ref. !, provides a framework in 
which this problem can, in principle, be addressed. The AdS/CFT 
correspondence relates a supergravity theory in the five- 
dimensional Anti-deSitter space AdS; to a strongly coupled, 
rank N, superconformal gauge theory on the 4-dimensional flat 
boundary R*' of AdS;. This is a holographic correspondence 
since R* is the (conformal) boundary of AdS;; hence it relates a 
gravitational theory in spacetime, to a conformal field theory 
(without gravity) on the boundary?~*. In particular, the classical 
geometry, ie., the supergravity vacuum, should therefore emerge 
from a particular quantum state of the dual gauge theory. The 
main aim of our work is to make this precise by exhibiting a 
canonical (i.e., background independent) such quantum state Vy 
and by showing that the supergravity vacuum in question 
emerges from the probability amplitude of in the t'Hooft limit 
where the rank N of the gauge group tends to infinity. 

In the general setting of minimal supersymmetry, the super- 
gravity vacuum is encoded by a Sasaki-Einstein metric gjy on a 
five-dimensional compact manifold M>°. On the gauge theory 
side the MN =1 superconformal symmetry is encoded by a 
complex cone Y of six real dimensions. This means that Y is a 
complex affine algebraic variety with a unique singular point yo 
(the tip of the cone), endowed with a repulsive holomorphic 
R, p-action, representing the conformal dilatation symmetry of 
the gauge theory. In the AdS/CFT correspondence the compact 
manifold M on the supergravity side arises as to the base of the 
complex cone Y of the gauge theory: 


M := (Y — {yo})/Rso. (1) 


The complex cone (Y, R,,) comes with a canonical holo- 
morphic three-form © which is homogeneous of degree 3. Such a 
space Y is often called a Calabi-Yau cone in the literature and is 
usually endowed with a conical Calabi-Yau metric gy, i.e., a Ricci 
flat conical Kahler metric, as well as a radial coordinate r, arising 
as the distance to the vertex point yo with respect to the 
Calabi-Yau metric gy. However, in the background independent 
formalism that we shall stress, the complex cone Y is not, a priori, 
endowed with any metric. This is crucial since our aim is to show 
how a metric emerges from the metric-independent BPS-sector of 
the gauge theory. Recall that BPS-states can be represented by 
holomorphic polynomial functions on the vacuum moduli space 
of the rank N gauge theory, whose mesonic branch is given by the 
symmetric product Y‘/Sy, where Sy denotes the symmetric group 
on N elements. The vacuum moduli space can thus be described 
in purely complex algebro-geometric terms’~?. 

In this work, we show that by imposing all the manifest 
symmetries of the rank N gauge theory one naturally arrives at a 
canonical quantum BPS-state Vy. More precisely, Yy can be 
realized as a wave function on the vacuum moduli space YN/Sy 
and its amplitude |‘¥|* induces a measure on YY: 


[Py(¥15  YyI?(Q.A Q)?™ (2) 


which is Sy-invariant and R.,,—invariant along each factor of 
YN (here Q denotes the conjugate of Q multiplied by —i so that 
QA Q is a volume form on Y). Somewhat surprisingly the state 
Wy does not appear to have been considered before in the lit- 
erature. Its explicit expression is given in Section II C, where the 
relation to the BPS-sector of the gauge theory is also explained. By 
quotienting out the conformal R.)—symmetry, we arrive at a 
measure on MN/Sy, which after normalization yields a canonical 
probability measure. Equivalently, we obtain a canonical ensem- 
ble of N “point-particles” on M. We show that a Sasaki-Einstein 


t is expected that a quantum theory of gravity should be able 


metric gy, on M emerges from the canonical ensemble in the large 
N-limit. In a little more detail, our construction goes as follows. 
First, the volume form dV, of gz emerges, after which the metric 
gu can be recovered from dV, in a standard manner (by simply 
differentiating dV, twice). In fact, in the course of this procedure 
an N-dependent radial function ry on Y naturally appears in an 
intermediate step and it induces a “quantum correction” gis) to 
the Sasaki-Einstein metric gj, on M. In this way, a limiting radial 
coordinate r on Y emerges from the gauge theory as N— o. 
From the perspective of the AdS/CFT correspondence, the radial 
coordinate r on Y corresponds to the radial coordinate on AdS; 
and thus our work reveals how the geometry of AdS-,, transversal 
to the conformal boundary R*!, naturally emerges from the 
gauge theory. A “spin-off effect” of this procedure is that it also 
produces (by a kind of back-reaction) a conical Calabi-Yau 
metric on Y, namely the cone over gy, whose distance to the 
vertex point yo is precisely r. It should be stressed that there are 
very few cases known where the Sasaki-Einstein metric gj, on M 
can be explicitly computed (but see ref. !° for a notable family of 
exceptions). Thus an important feature of our construction is that 
it furnishes canonical approximations gis) of Sasaki-Einstein 
metrics, encoded in terms of algebro-geometric data through an 
explicit integral formula. These can be numerically computed 
using Monte-Carlo methods, as detailed in section “Specialization 
to the toric case” for the toric case. Some relations to previous 
results and ideas for future work are discussed in section 
“Discussion”. 


Results 

Background: AdS/CFT and BPS-states. Recall that the low- 
energy dynamics of a general supersymmetric gauge theory is 
controlled by the moduli space of classical vacua M. The space 
M may be defined as the critical points, modulo complex gauge 
equivalence, of the superpotential W appearing in the (UV) 
Lagrangian of the gauge theory!!!*. Moreover, the BPS-operators 
of the gauge theory, i.e., the local operators preserving half of the 
supercharges, may be represented by holomorphic, polynomial 
functions on M. They thus form a ring known as the chiral ring 
of the gauge theory, which is graded by the R-charge. In the 
superconformal case, the BPS-operators can be viewed as states 
by the usual operator-state correspondence in CFT. The BPS- 
states are chiral primary states and saturate the BPS-bound, 


3 
A==R 3 
58 (3) 


where A denotes the conformal dimension. As a consequence, the 
BPS-sector tends to be robust under non-perturbative corrections 
and can thus be used to probe the strong-coupling regime of the 
gauge theory. 

In the setting of the AdS/CFT correspondence, the mesonic 
branch of the moduli space of vacua of the rank N gauge theory is 
the symmetric product® 


My = Syn" ¥ = ¥"/S,, (4) 


where Y is a complex cone. From the string theory perspective, 
this space parametrizes the transverse positions of ND3-branes 
inside Y. For example, in the maximally supersymmetric SU(N)- 
case, the superpotential W(Z), Z2,Z3) is defined on 3 complex 
NxWN matrices Z; (transforming in the adjoint representation) 
and W = Tr(Z,[Z,,Z;]). Thus the mesonic vacuum moduli 
space My may be parametrized by the set (cy /Sy of joint 
eigenvalues of (Z,,Z,,Z3),showing that, indeed, Y = C? in 
this case. 


The mesonic BPS-sector in the maximally supersymmetric case 
is isomorphic to the ring of holomorphic, polynomial functions 
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on (cy 5... Gauge theories with minimal supersymmetries 
may be constructed using quivers, encoding the matter content of 
the Lagrangian, as well as the superpotential!>. In general, this is a 
highly non-trivial task, but our approach only requires that the 
corresponding moduli space of classical vacua, encoded by the 
complex cone Y, is given. 


Mathematical prerequisites: complex-geometric setup. In this 
section we provide complex-geometric background, emphasizing 
a background-independent (i.e., metric-independent) perspective; 
see ref. '4 and the monograph!> for the more standard metric- 
dependent point of view. Let Y be a three-dimensional complex 
algebraic affine variety with an isolated singularity yo. Concretely, 
Y may be realized as the zero-locus of a collection of holomorphic 
polynomials on some complex space C™. Denote by J the 
induced complex structure on the regular locus Y — {yo}. Assume 
that Y is endowed with a 


e A repulsive holomorphic R , )-action that fixes yo 

e a holomorphic top-form Q, defined on the non-singular 
locus Y— {yo}, which is homogeneous of degree 3 with 
respect to the R.,)-action on Y. 


Such a space Y will here simply be called a complex cone, 
though it is often called a Calabi-Yau cone in the physics 
literature and an affine Gorenstein cone in the mathematics 
literature. The form © is uniquely determined up to a 
multiplicative constant and can often be written down explicitly. 

The vector field on Y generating the R., )-action will be called 
the dilatation vector field and denoted by 6. Rotating 6 with the 
complex structure J on Y yields another vector field that we shall 
denote by €: 


f= Jo. (5) 


The space O(Y) of holomorphic functions on Y decomposes with 
respect to the R,, y-action: 


BD O1.%), 


k=0,1,... 


where OQ, (Y) is the vector space of holomorphic (polynomial) 
functions, which are homogeneous of degree A; with respect to 
R,. From the perspective of the underlying superconformal 
gauge theory, the infinitesimal R.,-action 6 represents the 
conformal dilatation and 2& represents the R-symmetry (the 
factor 2/3 ensures that © has the same R-charge, 2, as the chiral 
superspace volume form d?6, where @ denotes fermionic coordi- 
nates of positive chirality). Thus, Eq. (5) is the complex-geometric 
realization of the BPS-relation (3) and the summands ©), (Y) in 
the decomposition (6) are thus BPS-states of dimension and R- 
charge equal to A,. 

The base of a complex cone Y is the compact five-dimensional 
manifold M defined by 

M = (Y— {y})/Ryo. (7) 

Thus, M is the base of the fibration (Y — {yy}) — M obtained by 
quotienting out the R., »—action on Y. Since the vector field € on 
Y commutes with the generator of the R.,,—action it induces a 
vector field on M that we denote by the same symbol é—known 
as the Reeb vector field on M in the mathematics literature. A 
metric gyz on M is said to define a Sasaki-Einstein metric on 
(M, €) if gy; has constant Ricci curvature, normalized so that it 
coincides with the Ricci curvature of the standard round metric 
on the unit-sphere, and is compatible, in a certain sense, with the 
complex structure on Y. 


The compatibility in question can be formulated in various 
ways, but the crucial point for our purposes is that a 
Sasaki-Einstein metric gj¢ can be explicitly recovered from its 
volume form dV, as follows. First observe that dV, induces a 
radial function, ie., positive function r on Y which is one- 
homogeneous with respect to the R , )-action: 


pany” 
_ ls A 8 
p= (208) ‘ 


where we have identified the volume form dv, gy OD M with its 
pull-back to Y. The metric gi on M with volume form dV, 


Sasaki-Einstein metric iff the corresponding radial fienction. ron 
Y solves the following PDE, after perhaps rescaling 1, 


(dd‘(r’))* =QA0, d :=J*d (9) 


on the regular locus Y— {yo} of Y. Here, d denotes the exterior 
derivative and d° denotes its “rotation” by the complex structure 
J, so that dd¢(r*) defines a two-form on Y and thus the three-fold 
exterior product (da‘(r?))° defines a six-form on Y. 

The PDE (9) is the celebrated Calabi-Yau equation on the 
complex cone Y;it is equivalent to the condition that the conical 
Kahler metric 


, isa 


gy = dd(ry(-,J), (10) 


on Y is a Calabi-Yau metric, ie, the Ricci curvature of gy 
vanishes. Finally, gy, may be explicitly recovered by identifying M 
with the level set {r= 1} in Y and letting gay be the restriction of 
gy to the level set {r= 1}. 

It is important to emphasize that, in general, the base (M, €) of 
a complex cone (Y,R,,) may not admit a Sasaki-Einstein 
metric?!4!6, Equivalently, this means that a complex cone may 
not admit a conical Calabi-Yau metric and hence no radial 
solution r to the Calabi-Yau equation (9). Indeed, by!® there 
exists a Sasaki-Einstein metric on (M, ) iff the complex cone 
(Y,R,,) is K-stable. This is a purely algebro-geometric 
condition. 

As shown in ref. !” the K-stability of Y can be viewed as a 
generalized form of the maximization condition for the a-central 
charge of the SCFT. This means that, in general, there are 
obstructions to the existence of a SCFT with a given IR , ,—graded 
mesonic chiral ring O(Y“/Sy). In the present approach, a 
different, but conjecturally equivalent, stability type condition 
naturally appears, which is a variant of the notion of Gibbs 
stability introduced in the context of Fano manifolds in ref. 1°. 


Main results: emergent geometry. According to the AdS/CFT- 
correspondence the classical supergravity vacuum geometry in 
AdS; should emerge from some limit of a quantum state in the 
dual CFT gauge theory on the boundary. Concretely, in the 
present setting the non-trivial part of the supergravity vacuum in 
question is encoded by a Sasaki-Einstein metric gy on the 
internal compact space M, corresponding to the base of the 
complex cone Y of the dual gauge theory®. Hence, we want to 
show that the Sasaki-Einstein metric gy on M emerges in a 
certain “large N-limit” of a specific (and background free) 
quantum BPS-state ‘Yy in the dual-rank N gauge theory. 

First, recall that Y is endowed with a holomorphic top-form Q 
and hence one can endow the mesonic classical vacuum moduli 


space YN/Sy with the volume form (QO. A 0)”. Let Wis Vy, 
be a maximal number of linearly independent mesonic BPS-states 


for the rank 1 gauge theory with the same R-charge A,. In other 
words, Y\,..-,Wy, form a basis in the space O,(Y) of 


holomorphic functions on Y with the same charge A,. Denote 


NATURE COMMUNICATIONS | (2022)13:365 | https://doi.org/10.1038/s41467-021-27951-9 | www.nature.com/naturecommunications 3 


ARTICLE 


NATURE COMMUNICATIONS | https://doi.org/10.1038/s41467-021-27951-9 


by Wace the corresponding Slater determinant, ie., the totally 
antisymmetric holomorphic function on Y'« given by 


Fdet Wr + n,) = 1, (1) YeayV) Woy On,): (11) 


The function Y ge is independent of the choice of bases in O71,(¥) 
up to an overall multiplicative constant. It thus defines a baryonic 
BPS-state in the rank N; gauge theory®!?. 

We aim to construct a symmetric measure on the mesonic 
classical vacuum moduli space Y%* / Sy,» Which is invariant under 
the conformal R,,—action on each factor. One might be 
tempted to try with the density |Yge|*, but this is unfortunately 
not (IR ,,)*-invariant. The resolution is to simply take a suitable 
fractional power of the Slater determinant, namely 


—3/A, 
= We (12) 


Then, an (IR,,)‘*-invariant measure on YN/Sy is given by the 
combination 

[iy One ay) AQ) (13) 
This measure is both invariant under the Sy, permutation 
symmetry and the conformal R-symmetry, as well as the 
R-symmetry. However, since Y is non-compact, the integral of 
this measure over YN« / Sy, diverges. To circumvent this problem, 
we simply quotient by the (R.,)‘*-action to get an induced 
measure on the compact space M« /Sy,- To be specific, 
contracting the top-form QAQ on Y with the dilatation 
vector field 6 yields a 5-form 1;(QAQ) on Y. Thus, the 
(IR, ))**—invariant form 


en Oia? = Yn JI sQ.A O) (14) 


maybe identified with a measure on the compact quotient space 
MN: /Sy,. Finally, the canonical probability measure dP, on 


MN« / Sy, is defined by 


@N, 
’ 


1 —_ 
dP = Zz [Wy, I? (49(O A Q)) (15) 
k 


where 


ae (16) 


Zn = Wy, 1? (4g(Q A Q) 
k Me cn k 

Note that in (15) we have made the implicit assumption that 
Zy, is finite. Since |Vy, |? blows-up along a hypersurface in 
YN« / Sy,» namely the zero-locus of the Slater determinant 
Waep this is actually a very non-trivial condition. We interpret it 
as a consistency condition (which turns out to be related to the 
mathematical notions of K-stability and Gibbs stability). 

Since Vy is holomorphic away from its singularity locus in YN/ 
Sw (the vanishing locus of ‘¥4,,) the state Vy can be viewed as a 
bound state of BPS-states in the rank N gauge theory. 

We emphasize that by “canonical” we mean that the definition 
of dPy, is background independent, in the sense that it does not 
depend on any underlying metric on Y or M. It only depends on 
the complex structure J on the classical vacuum moduli space and 
the R.)-action and thus only on the superconformal symmetry 
of the rank N-gauge theory. This is a crucial point for describing 
the emergence of the classical Sasaki-Einstein metric gy on M, 
which is our main focus. 

We have also restricted the values of the rank N to be a 
sequence of integers 


N, := dim, ,(¥), (17) 


ie, the multiplicity of the R-charge A,.This can be seen as a 
quantization condition. As is well-known, in the quasi-regular 
case (discussed in the next section) N; is a polynomial in k of the 
form, 
2 

N= ZV + Ol), h4~k, k= co (18) 
where the positive number V is an algebraic invariant of the 
complex cone (Y, R,.), known as its volume?!*+!6, 

Assume, for simplicity, that the complex cone (Y,R, 9) 
associated to the gauge theory is quasi-regular. A complex cone 
(Y, R,,) is quasi-regular if (up to a rescaling) the R.,)—action 
on Y can be complexified to a holomorphic C* —action. Denote 
by dPY the probability measure on M defined as the 1 — point 
correlation measure of the canonical ensemble (dP, M% /Sy) 


introduced in the previous section. In other words, dP? is 
obtained by “integrating out” all but one of the factors of MN: 


1 
ary) == | pense, EMO Is oo IOP 


N 


(19) 
x (1,(0.A D)) PN". 
Our main conjecture can now be stated as follows: 
Conjecture A: Assume that the canonical ensemble (dP, MY / 
Sy) is well-defined, ie., that Zy <oo. Then 


(i) The one-point correlation measure dP? converges, as 
Nc, to the volume form dV, of a Sasaki-Einstein 
metric gy on (M, é), normalized to have unit-volume, 


Jim dP = dVx; (20) 
(ii) The sequence of radial functions 
a) ~1/6 
a (ee ae (21) 
N 
13(0 A QO) 


on the complex cone (Y, R,9) converges, as N — oo, toward 
the radial function r of the Calabi-Yau metric gy on Y 
corresponding to the Sasaki-Einstein metric gy, on (M, ). 
Conversely, if there exists a unique Sasaki-Einstein metric 
gu on M, then Zy < oo. 


(iii) 


We emphasize that for finite N, the radial function ry yields an 
explicit approximation oe to the Sasaki-Einstein metric gj on 
M by identifying M with the level set {r,;;= 1} and setting 


gyn = dd (ry) CsI) 


(c.f. Eq. (10)). Hence, part (ii) of the conjecture is equivalent to 


(22) 


(N) 


Jim gu’ = &™ 


(23) 


which is the sought-after emergence of the Sasaki-Einstein metric 
on M in the large N limit. 

To be mathematically precise, the convergence statements in 
the conjecture are supposed to hold in the standard weak 
topologies. In fact, we make the stronger conjecture that the 
random measure N7! as 6, on the canonical ensemble 
converges in law toward the deterministic measure dVjy. Below 
we will prove a 6-deformed version of this conjecture. 

We first introduce a real-analytic family of probability 
measures dP g on M\N/Sy, defined for a real parameter , such 
that dPy g coincides with dPy for B= —1, if Zy <ov. To this 
end, fix a background radial function rg on Y. We can then 
identify the base M := (Y — {yo})/IR,9 of the cone Y with the 
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level set {re = 1} in Y and define dP, as follows: 


__} 2g y@N 

dPyg = aan Wace 1s Vo» Yn) kdVo ’ (24) 
NB 

where dVo denotes the volume form on M obtained by restricting 
the five-form 1;Q A Q to the level set {% = 1} and Zy_p is the 
corresponding normalization constant (recall that N is the 
multiplicity of the charge A;). The parameter 6 can be viewed 
as a regularization parameter, since Zyg is automatically finite 
when $ > 0 (or slightly negative). However, it should be stressed 
that it is only in the canonical case 6 = —1 that the probability 
measure (24) is independent of the choice of radial function 
ro. Let ry.g be the radial function on Y defined by 


(25) 


(coinciding with rj when B=—1) and denote by Bios the 
corresponding metric on M, obtained by replacing the radial 
function ry in Eq. (22) with ry,g. We then have: 


Theorem B: For each f > 0, there exists 
(i) a volume form yg on M such that 
: GQ) _ 
i AE Ns = Hp 
(ii) a radial function rg on Y such that 


lim ry,=r 
N->0oo NB B 


(iii) and a Sasaki metric guyz on M such that 
lim gO), = gy 9. 
Noo MB NB 


Moreover, if (M, €) admits a Sasaki-Einstein metric, then rg and 
8g extend real-analytically to [—1, 9] and setting B = — 1 yields 
a Sasaki-Einstein metric gy, on M. 

The proof is given in SIV. In the course of the proof, we will 
show that the square of the limiting radial function rg is the 
unique conical Kahler potential on Y solving the following PDE 


on Y — {yo}: 
2 3(B+1) 
(dd°r3)’ = (") QAQ. 
1G 


0 


(26) 


In particular, for B= -—1 this is indeed the Calabi-Yau 
equation (9) for the radial function r corresponding to a 
Sasaki-Einstein metric gj, on M. 


Finiteness properties of the normalizing constant. Loosely 
speaking, Theorem B thus shows that Conjecture A holds after 
the analytic continuation. More precisely, it shows that Con- 
jecture A holds under the assumption that the order of taking the 
limits N,— co and B — —1 may be interchanged. By a physics 
level of rigor Conjecture A may thus be considered as established. 
However, we do expect that the introduction of the parameter f is 
not needed and, in particular, that Zy < oo if and only if (M, 4 
admits a unique Sasaki-Einstein metric. In the case when the 
Sasaki-Einstein metric is not unique, ie., when the Lie algebra 
a(Y, &) of the automorphism group of (Y, is non-trivial?°, we 
conjecture that Zy ,<oo for any B>—1 when N is sufficiently 
large. The “only if direction” can be deduced from recent 
mathematical results in complex geometry for complex cones Y of 
any dimension, as will be shown elsewhere. Proving the 
remaining direction appears, however, to be very challenging, 
except in the case when Y has complex dimension two, where a 
direct proof of Conjecture A can be given. 


For example, when Y = C? realizing M as the Hopf fibration 
over the two-sphere S? and factorizing Vj,(x,, x5, .--; Xj) reveals 
that Zyg can be expressed as the configurational partition 
function for N particles on S? interacting by the 2D-gravitational 
force with a mean-field scaling: 


2B 
= _ x 1g Aen 
Zyp = Cy ie . II lx; — xlRsdA°™, 


(27) 
) 1<i#j<N 

expressed in terms of the restriction to S* of the Euclidean norm 
on R° (cf. Eq. (36)). Applying the arithmetic-geometric means 
inequality reveals that the integral is finite precisely when 
B>-—(1—1/N).A similar argument applies to any Y of complex 
dimension two, using that Y is a Kleinian singularity, ie, Y= 
C?/G for a finite subgroup G of SU(2) and thus that M is a Seifert 
fibration over S?,branched over three points (in this case 
Zy.p < co for B =—1-when N is sufficiently large-since 9(Y, ¢) 
is trivial; details will appear elsewhere). For higher-dimensional Y 
the Slater determinant V4.,(x,,X,..., xy) can not, however, be 
factorized. But a condition ensuring that our canonical partition 
function Zy is finite for any N is that Y is an exceptional 
singularity?2. For example, there are exactly (up to conjugation) 
five cases of exceptional (non-isolated) singularities of the form 
Y = C°/G, for G a finite subgroup of SL(C, 3); notably Klein’s 
simple group of order 168, PSL(2,7),”!, Cor 3.15. See ref. 73 for 
the construction of the matter content and gauge groups of the 
corresponding gauge theories; the quiver graphs for the five 
“exceptional” groups G appear in figure 5 in refs. 7°, 

Moreover, a list of three-dimensional exceptional quasi- 
homogeneous (isolated) hypersurface singularities in C* is given 
in ref. 2*, Cor 1.1. Consider for example the case when Y is a 
Briskorn—Pham singularity: 

{Ze + Zp +Z74+ZP =0} CC, (28) 
endowed with the diagonal C*—action with weights proportional 
to a; '. If the powers a; are coprime and taken in increasing order, 
then Y is an exceptional Gorenstein affine variety iff 
1< ,a;!<1+4a;!. This condition is, for example, satisfied 
for the powers (2,3,7,11) and is, in fact, equivalent to the 
condition that Z, < oo for any N. It also implies that the base of 
Y admits a Sasaki-Einstein metric. However, in Conjecture A we 
only demand that Zy <oo for N sufficiently large. 


Discussion 

Let us conclude by briefly mentioning some relations to previous 
work. First of all, our work is very much in the spirit of the 
program for emergent geometry in AdS/CFT initiated by 
Berenstein’+*°. The main new feature in our work is the 
appearance of a negative and fractional power of the Slater 
determinant in the definition of the state Vy (see Eq. (12)) and its 
B-deformation. This is crucial in order to obtain background 
independence and to see the emergence of the spacetime metric, 
as explained in Section II C. 

As explained at the end of section “Methods”, our approach 
builds on the probabilistic approach to Kahler—Einstein metrics 
on Fano manifolds introduced in refs. !87°7, which, in turn, is 
motivated by the Yau-Tian-Donaldson (YTD) conjecture for 
Fano manifolds. A different connection between the YTD con- 
jecture and AdS/CFT was first exhibited in ref. 17 (compare the 
discussion on stability in section “Mathematical prerequisites: 
complex geometric setup”). 

Finally, we note that our canonical ensemble on M may be 
viewed as an ensemble of N dual giant gravitons, as will be 
elaborated on in a separate publication. 
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Methods 

Proof of Theorem B. We will show how to deduce the theorem from the results in 
refs. *”2° concerning a probabilistic approach to Kahler-Einstein metrics on Fano 
manifolds. We thus start with some well-known geometric preparations to realize 
Mb, as a fibration over a Fano manifold (orbifold); see ref. 3°, Section 2.3. First 
assume that € is a regular Reeb vector field. This means that the orbits of the 
complexification of & coincide with the orbits of a C*—action on Y without fixed 
points on Y*: = Y — {yo}. The corresponding compact complex manifold X := 
Y*/C* is a Fano manifold, ie., the dual K% of its canonical line bundle Ky is 
ample. The natural projection from Y to X realizes Y* as the total space of the qth 
tensor power Ke! — X for some rational positive number q, when the zero-section 
has been removed (and Y gets identified with the variety obtained by blowing down 
the zero-section). For example, in the cases when Y is C? and the conifold one gets 
X = P* and X = P'x P' with q=1/3 and q=1/2, respectively. In other words, 
denoting by L the ample line bundle (K%)®4, we may identify Y with the total space 
of the fibration L* — X, with the zero-section deleted. The fixed radial function ro 
on Y corresponds to a Hermitian metric ||-|| on L* and the induced quotient 
fibration My — X realizes Mo is a principal U(1)-bundle over X, namely the unit- 
circle bundle of (L*, ||-||): 


M, © L*(=Y*) 


\ 1 (29) 
4 
where 
My = {rm =1} =f - =U. (30) 


As a consequence, there is a one-to-one correspondence between the space P(M,)* 
of é-invariant probability measures 4 on Mp and the space P(X) of probability 
measure v on X: 


(31) 


expressing yu as the fiber product of v with the é-invariant probability measures dO 
defined on the fibers of the fibration My — X. In other words, v is proportional to 
the contraction of y with & descended to X. Introducing local holomorphic coor- 
dinates z on X and locally trivializing L* with the holomorphic section (dz)®4 of 


Ke! we may locally express 


u=vedé, 


Q= dea dw'/9), n= (|w|2et@) 794, (32) 


where w is a local holomorphic coordinate along the fibers of L* and e*) denotes 
the squared norm of (dz)®4, i.e. eo) =|| (dz)®4||?. The local formula for Q fol- 
lows from the observation that dz d(w!/2) glues to define an equivariant global 
holomorphic three-form on Y*. The appearance of the power 1/3q in the formula 
for r2 then follows from the relation & = 3qé,., where is the standard U(1)- 
action along the fibers of L* (satisfying €,.w = iw), resulting from the normal- 
ization condition ¢Q = 30. Since the weight space H),(Y) may be identified with 
the space H® (X, L®*) of holomorphic sections of the kth tensor power of the 
holomorphic line bundle L — X it also follows that A; = 3qk. Concretely, the 
identification in question is obtained by noting that an element ¥ in H, (Y) may be 
locally expressed as ¥(z, w) = filz)wk for a local holomorphic function f,(z), 
globally transforming as a holomorphic section of L®* — X. In particular, 


Pact =f aet(Z1 + --5 2N JL Wy (33) 


where the local holomorphic function f j,,(Z,, -.-;Zy) on X% transforms as a 
IN 


holomorphic section of (L®*)"" — X%, namely as the Slater determinant of 
H(X, L®*), Moreover, by Eq. (32) we have 


: 2ekeule) ... e-kfalew) 
| act ney Uf det(Z1; sey Sy) [OPO oe ko (Zn és 
= I faell?, 
using, in the last equality, the induced Hermitian metric ||-|| on (L®*) Ne _, XN 
With these preparations in place, we can thus express 
dP y's = yg @ d0°%, (35) 
where vy, is the probability measure on XN defined by 
28 
kad @N 
ve aM faalltaV ee 


Np — 28 5) 
Sxe l| Factll*#dV xe 


where the volume form dV on X corresponds to the volume form dVp on Mo 
under the correspondence (31). The probability measure vy,g on XN is precisely the 
probability measure defined by the “temperature deformed” determinantal point 
process on X introduced in ref. 7’, associated with the Hermitian holomorphic line 
bundle (L, ||-||) over the compact complex manifold X endowed with a volume 
form dVy at the inverse temperature $/q. By ref. 2”, Theorem 5.7 (and ref. 27, 


Lemma 5.1) its one-point correlations measures We converge as N— ©, in the 
weak topology of measures on X, toward a volume form vg on X of the form 


B 
Vp= e** dV x for the unique smooth function gg on X satisfying the following PDE 
on X: 


= (w + dd'g,) =e dV x, (37) 


where w,, is the Kahler form defined by the curvature of the metric ||-|] on L, locally 
expressed as w, = dd°$o(z). Using that w, = dd‘log (rg) /3q a direct calculation 
now reveals that the radial function rg on Y defined by 


1/6, 
ee fo a ade 
Be \dVo - 


satisfies the PDE (26). This proves the convergence in item 1 of Theorem B with yg 


(38) 


the volume form in P(M)é corresponding to vg in P(X). Similarly, applying?’, 
Corollary 5.8 then proves the convergence of ry,g and gy,g toward rg and gp, 
respectively. Moreover, if (M, &) admits a Sasaki-Einstein metric, then it corre- 
sponds to a Kahler—Einstein metric wy on X. Thus, as shown in Step 2 of the proof 
of ref. 2°, Theorem 7.9, the PDE (37) on X admits a unique solution gp for any 
B>—1 and 9g defines a real-analytic family converging to the Kahler potential g_, 
of a Kahler-Einstein metric on X,as 8 — —1. When rephrased in terms of rg and 
&ug this concludes the proof of Theorem B in the regular case. Finally, in the case 
when € is quasi-regular one can proceed essentially as in the regular case, using that 
in this case, the quotient X is a Fano orbifold, so that the role of Ky is now played 
by the orbifold canonical line bundle of X. 

This concludes the proof of Theorem B. 


Specialization to the toric case. In this section, we specialize our proposal to the 
case of a toric quiver gauge theory. As is well-known in this case the corresponding 
complex cone Y is a toric affine Gorenstein variety. As shown in refs. >! (Y, &) 
admits a conical Calabi-Yau metric gy iff is the unique minimizer of the volume 
functional V(é) on the space of normalized Reeb vectors, introduced in ref. 39. 
Equivalently, from the gauge theory point of view, this means that the U(1)r- 
symmetry induced by € maximizes the a-central charge. We explain how gy 
emerges from our proposal using a tropicalization procedure, which renders the 
proposal provably convergent and computationally feasible. It also applies to 
irregular Reeb vector fields. 

Let Y be a 3-dimensional normal affine toric variety. This means that Y is a 
normal affine variety endowed with the holomorphic action of the 3-dimensional 
complex torus 3 with an open dense orbit, where T;, := (C*)* denotes the 
complexification of the compact torus T: = U(1)°. Accordingly, we can identify T, 
with an open subset of Y and view Y as an T;—equivariant compactification of 
Tc. Denote by yo the unique point in Y which is fixed under T; and by 
Z= (2, 22, Z3) the holomorphic coordinates on T.. Since it requires no extra effort 
we will allow yo to be a non-isolated singularity, following*. 

The ring R(Y), consisting of all holomorphic polynomials on Y, splits with 
respect to the T;-action on Y: 


RY)= @ Ce, 


pel*nZ 


a= aigzy, (39) 


where C* is the “moment polytope” of the affine toric variety Y. C* can be 
represented as the convex cone in (R°)* whose dual is a convex cone C Cc R?. To 
simplify the notation we will identify the dual (IR*)" with R° in the usual way. The 
Reeb vector fields £ on Y may be identified with vectors & in IR? lying in the interior 
of the cone C. Denote by A, the corresponding weights (&,p) as p ranges over 
C* Z?, ordered so that Ay <Ap...<.We can thus represent the corresponding 
weight spaces of R(Y) as 

Nx 


H,(Y) =@ Ce, 


i=1 


(40) 


where we have enumerated the lattice points p,,..., Py, in the 2-dimensional 
convex polytope P; defined as the intersection of the convex cone C* with the 


hyperplane {(&,-) = A,}: 


Prh=APpOZ?, Pee {CN (E,-) =}, (41) 


(note that the polytope Pr is denoted by A in ref. 9 and called the characteristic 
polytope). Now assume that Y is Gorenstein and denote by © the T,—equivariant 
holomorphic top-form on Y — {y°}.On Tx C Y, we can express 


tg BE Oe 


aQ=2'0,, 
Z 22 23 


(42) 
for some I € Z, where Qp is the standard T,—invariant volume form on Tc. The 
condition that Q is homogeneous of degree 3 under the Reeb field & translates into 
the condition (I, ) = 3. For example, when Y = C3 with the standard Reeb vector 
& = (1,1, 1), the cone C* is the positive octant in R? and P; bounds the unit- 
simplex. We will also briefly discuss the case of the conifold below. 

We now specialize our proposal to the toric case. First, note that using the 
bases z?1,...,2?% in H. y,(¥) the corresponding Slater determinant '¥., may be 
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represented as the following holomorphic function on Te c YN: 


Y ae = — lel Pou... zPony) 
det(Z1 5 12y,) a yz Poin, (43) 


Accordingly, on the open dense subset T’ of YN we can, using formula (42), 
express 


@N, Nx 


Py, PQA OY = py (B55 ZW AO)” * (44) 


where 


23. 
Ak 


>» (=1)!4! zea) wee Zloty) ; (45) 


FESy,, 


Py, Ziy Zy,) = 


Here we have defined q;: = pi/A, — 1/3, i= 1, ..., Nj, corresponding to the discrete 
points of the scaled and shifted polytope Q, defined as 


Qi = Py/Ay — 1/3. C RP N(Z/A,)’ 
Qe = P, — £/3, 


where Q; is the limit of Q, when k — o». 
From a computational point of view, this explicit expression for py (2), .-.,2y,) 


(46) 


is still rather challenging to work with directly. But the construction of a 
Sasaki-Einstein metric can be simplified by further leveraging the toric structure. 
To see this, first recall that, in general, the group G(Y, €) of all biholomorphisms 
of a complex cone Y, commuting with the Reeb vector field § and homotopic to the 
identity, acts transitively on the space of conical Calabi-Yau metrics gy”. In 
particular, the toric case G(Y,&) contains the group T. and thus gy is not uniquely 
determined, but can be taken to be T-invariant. The density py, (2, ...,y,), on the 
other hand, is not TN:—invariant. This is to be expected as the large N-limit should 
encapsulate all conical Calabi-Yau metrics on Y—not only the T-invariant ones. In 
order to directly extract a T-invariant conical Calabi-Yau metric gy from the large 
N-limit we can modify the density py, on Y%« to render it TN*—invariant. This can 
be achieved in various ways, but from a computational point of view, the most 
efficient modification appears to replace the density py, (2, ..-,Zy,) with its 


8 
2 
) ; (47) 


where, in particular, the sum over Sy, has been replaced with a maximum. In terms 
of the logarithmic real coordinates 


X = (1, X),%3) = (log (|z, 17), log (z,|7), log (|z3|7)) eR (48) 


this means that 


tropicalization: 


gin ne ain.) 


(Nx), : 
Pirop (21 + Zy,) = | Max 
oeSy, 


(N,) — NENW (x, ..., 1) 

Prrop (215-2, ) =e BNE (49) 

where ie (x1, .-.,Xy,) denotes the following symmetric piece-wise affine convex 
P N 

function on (IR?) *: 


3 
(N;) = 
E trop (X13 +1 Xy,) = Ne (1-400) pet 


an (25,:4o0)))- 


Hence, the corresponding T*—invariant “tropicalized” measure on YN* may be 
expressed as follows 


(50) 


pov (yy ns 2y, MOA OM = Neier )(deJOME @ (AB)EX*, (51) 
with d@ denoting the T-invariant probability measure on T. Contracting Q A O 
with the dilation vector field 5: = —Jé, as before, thus yields a TN: invariant 
measure on Ms. By performing a linear change of coordinates on R°, we may 
assume that dilatation on Y corresponds to translations in the x3-variable in R*. 
The contraction in question thus corresponds to replacing R? in Eq. (51) by R’. 
However, the integral over M%« given by the corresponding normalizing constant 
always diverges due to the diagonal action of the residual symmetry group T,./T. 
This action has the effect of translating the center of mass in IR? of a configuration 
221 
(x1, .-,Xy,) € CR“) 
be achieved by introducing a background radial function (as in Theorem B). But 
from a computational point of view the most efficient way is to simply impose the 
constraint that the center of mass in R? of (x,, +) Xy,) vanishes, ie., that x, + 
et Xy, vanishes. Finally, since the number N; may scale as o(k2), unless & is 


‘The remedy is to break the symmetry in question. This can 


quasi-regular, we replace N; with any positive integer N and take qj, ..., qu to be any 
sequence of points in the polytope Q; with the property that 


1 N 
lim — > 54, = VQ (52) 


Noo N j=1 


where Vg, is the Euclidean measure on Qe normalized to have unit total mass. In 


summary, the tropicalized probability measure thus corresponds to the following 
Boltzmann-Gibbs measure on (R?)” : 


1 (N) 
(N) —NE™) (20), Hi @N 
Htrop “= Tony & oe *w)(dx) , (53) 
trop 
with 
(N) 
Ze) (eo ak Cee ¥y) dy @N (54) 


together with the constraint of vanishing center of mass. 

We then have the following result. 

Theorem C: Ze. <oo for N sufficiently large iff (M, &) admits a conical 
Sasaki-Einstein metric (i.e., iff € minimizes the volume functional V(&) on the Reeb 
cone). Moreover, if this is the case then the law of the empirical measure capil 1 5x, 
on the ensemble (eR, 1, converges in law, as N — , toward a volume form 
Hse on R?; the normalized volume form of the unique T-invariant Sasaki-Einstein 
metric on (M, é) with vanishing center of mass (when expressed in real logarithmic 
coordinates on R?). 

The probability measure js; on IR? may be expressed as 


Usp = @ *Pdx, (55) 


where the function ¢ on IR? corresponds to the Kahler potential r? of the 
corresponding T-invariant conical Calabi-Yau metric on Y, ie. 7 = e?+I)/3, 
Theorem C thus provides evidence for Conjecture A. The advantage of the 
tropicalized setup is that the corresponding energy Euep(X 15 ..-,Xy) is continuous. 
The theorem can be shown using results in ref. 34, where a closely related tropical 
approach to Kahler-Einstein metrics on toric Fano varieties was introduced. 
Details will appear elsewhere. 

For example, in the homogeneous cases of C? and the conifold it is well-known 
that the polytope Q; is a translation of the two-dimensional simplex or the unit- 
square, respectively. In these cases, the function ¢(x) is given explicitly, modulo an 
additive constant, by ¢(x) = log(1 + e* +e) — x,/3 — x,/3 for C® and ¢(x) = 
log (e~*:/? + e*:/?) 4 log (e2/? + e*2/) for the conifold. 

Similarly, in the case when Y is C” the polytope Q¢ is equal to [—1/2, 1/2] and 
the corresponding one-point correlations can be computed explicitly for any finite 
N. The result is a polynomial of degree N in e~|, converging, as N — ©, toward e 
~29(), where $(x) = log(e~*/? + e*/?) — C. Indeed, in this case, the corresponding 
tropical energy fee turns out to coincide with the energy of a self-gravitating 
system in 1D (with a mean-field scaling), to which the exact results in ref. 3° apply. 

In general, however, the solution ¢ can not be explicitly computed. An 
important feature of Theorem C is that it provides an efficient way of obtaining 
numerical approximations to the solution ¢ and thus to the Sasaki-Einstein metric 
on M, using Hamiltonian Monte-Carlo. The starting point is the observation that 


the gradient of the energy NEW appearing in the Boltzmann-Gibbs measure (51) 
is, for a generic configuration (x,,...,x,) € IR’, precisely the discrete optimal 
transport map matching the points (x, ...,x,) with the fixed points qj, ..... qn in 
Q:*°. Thanks to the last years rapid developments of scalable optimal transport 
solvers this allows one to numerically compute the gradient in nearly O(N) 
operations (efficiently implemented on GPU-hardware*”3*), Moreover, since EW, 
is convex Hamiltonian Monte-Carlo should merely require O(N"/*) gradient 
evaluations*®, suggesting that the total running time of the simulation nearly scales 
as O(N!+1/4), 

Monte-Carlo techniques have previously been applied in ref. 4° to the vacuum 
moduli space YN/Sy in the case Y = C?, but using a different N-particle BPS-wave 
function (with a Metropolis algorithm). 

We leave the implementation of our simulation scheme for the future. 
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